{
  double PEnergy[45] = 
	{ 2.96171, 3.12387, 3.28378, 3.44595,	
	  3.60586, 3.76577, 3.92793, 4.09009,	
	  4.25225, 4.41441, 
	  4.57432, 4.73423, 4.8964 , 5.05631,	
	  5.21847, 5.37838, 5.53829, 5.6982 , 
	  5.85586, 6.01351,	
	  6.17568, 6.33559, 6.4955 , 6.65541,	
	  6.81532, 6.97523, 7.13514, 7.29505,	
	  7.45721, 7.61486,	
	  7.77252, 7.93243, 8.08784, 8.24324,	
	  8.4009 , 8.55405, 8.71396, 8.86937,	
	  9.01802, 9.17568,	
	  9.32883, 9.48649, 9.63964, 9.79054,	
	  9.9527  
	}; 

   double press = 101;//kPa
   for(int i=0; i<45; i++) 
   {
       double wave = 1240/PEnergy[i];
       double ref = 1 + (0.124523e-6/(2.61154e-4 - pow(wave, -2)))*press/101;
   
       cout<<i<<" "<<PEnergy[i]<<"*eV "<<ref<<endl;
   }
}
